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Energy transport in weakly nonlinear wave systems 
with narrow frequency band excitation 

Elena Kartashova* 
Institute for Analysis, J. Kepler University, Linz, Austria 

A novel discrete model (D-model) is presented describing nonlinear wave interactions in systems 
with small and moderate nonlinearity under narrow frequency band excitation. It integrates in a 
single theoretical frame two mechanisms of energy transport between modes, namely intermittency 
and energy cascade and gives conditions when which regime will take place. Conditions for the 
formation of a cascade, cascade direction, conditions for cascade termination, etc. are given and 
depend strongly on the choice of excitation parameters. The energy spectra of a cascade may be 
computed yielding discrete and continuous energy spectra. The model does not need statistical 
assumptions as all effects are derived from the interaction of distinct modes. In the example given 
- surface water waves with dispersion function lo 2 = g k and small nonlinearity - D-model predicts 
asymmetrical growth of side-bands for Benjamin-Feir instability while transition from discrete to 
continuous energy spectrum excitation parameters properly chosen yields the saturated Phillips' 
power spectrum ~ g 2 cu~ 5 . D-model can be applied to the experimental and theoretical study of 
numerous wave systems appearing in hydrodynamics, nonlinear optics, electrodynamics, plasma, 
convection theory, etc. 
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I. INTRODUCTION 

A central topic in the theory of weakly nonlinear wave 
interactions is the mechanism of energy transport be- 
tween modes. Considering what we can describe in the- 
ory and observe in experiment, there is good reason to 
believe that in any weakly nonlinear dispersive wave sys- 
tem there are two main types of energy transport: inter- 
mittency which is a periodic or chaotic exchange of en- 
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ergy among a small number of modes, and energy cascade 
which is a unidirectional flow of energy through scales in 
Fourier space. 

In systems with distributed initial state energy trans- 
port is studied in the frame of kinetic wave turbulence 
theory (WTT) by means of the wave kinetic equation, 
[1, 2]. In this paper we explicitly study wave systems 
with narrow frequency band excitation. 

Intermittency is based on finite size effects in a res- 
onator. The general properties of weakly nonlinear wave 
systems showing intermittency have first been character- 
ized through the solution of the kinematic resonance con- 
ditions, [3], which reflect the geometry of the resonator. 
The general dynamical characteristics of this type of en- 
ergy transport have been studied in the frame of discrete 
WTT, [4], for systems with narrow frequency band exci- 
tation. Main mathematical object of the discrete WTT 
is a set of dynamical systems for the amplitudes of in- 
teracting waves; each dynamical system corresponds to 
a resonance cluster composed of a small number of reso- 
nant triads or quartets having joint modes, [5]. 

Energy cascades in systems with narrow frequency 
band excitation have recently been described in [6] us- 
ing increment chain equation method (ICEM) . An energy 
cascade is represented as a chain of modes with nonlinear 
frequencies triggered by modulation instability at each 
cascade step. The energy spectra E(uj) obtained by the 
ICEM have exponential decay and can be written as: 

i>2 
i=l 

where for given linear dispersion function uj ~ k a , Ci are 
known functions of excitation parameters and 7$ vary for 
different magnitudes of nonlinearity. For comparison, in 
the systems with distributed initial state, studied in the 
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frame of kinetic wave turbulence theory (WTT), energy 
spectra decays according to a power law, 



E(lo) <~ lo 7 , 7 > 



(2) 



with different 7 for different wave systems, [1, 2]. 

In this paper we present, based on the resonance con- 
ditions, a common mathematical model, called D-modcl 
("D" for "discrete"), incorporating both forms of energy 
transport, intermittency and cascades, and give criteria 
under which conditions to expect which behavior. 

In D-model, intermittency occurs for very small non- 
linearity, < e < 0.1, provided that the geometrical form 
of the resonator permits resonance. An energy cascade 
occurs at larger levels of nonlinearity, e <~ 0.1 0.4, and 
its spectrum does not depend on shape or finiteness of the 
interaction domain. The outcome of the model strongly 
depends on the excitation parameters. 

D-model can explain the following phenomena ob- 
served in systems with narrow frequency band excitation: 

- no cascade but recurrent wave patterns are observed, 
[7] (surface water waves); 

- a cascade consisting of two distinct parts - discrete 
and continuous; form of spectra does not follow a power 
law, [8] (thin elastic steel plate); [9] (gravity-capillary 
waves in mercury); 

- a discrete energy cascade develops a strongly non- 
linear regime yielding breaking, a continuous part of the 
spectrum is not observed, e.g. [10] (surface water waves); 

- form of energy spectra depends on the parameters 
of excitation, e.g. [11, 12] (gravity surface and capillary 
water waves correspondingly); 

- amplitudes of direct and inverse cascades are not 
symmetric, e.g. [13-17]; 

- interactions of waves over several orders of magnitude 
[18] (capillary waves in helium). 

The model is briefly described in Sec. II. To demon- 
strate how D-model works, we give an example determin- 
ing cascade direction and scenarios of cascade termina- 
tion for surface water waves, depending on the excitation 
parameters (Sec. III). 

In Sec. IV we compare assumptions and predictions of 
D-model and kinetic WTT to give an experimentalist 
clues which model to apply in a given experimental set- 
up. A short list of conclusions and open questions is given 
in Sec.V. 



A. Here A, k, u> = u>(k) denote amplitude, wavevec- 
tor and dispersion function correspondingly. The small 
parameter is usually introduced as wave steepness e = 
Ak, k — |k|. If the nonlinearity is small enough, only 
resonant interactions have to be taken into account. The 
resonance conditions read 



for 3 waves: 



for 4 waves: 



w(ki) + w(k 2 ) = w(k 3 ), 
ki + k 2 = k 3 . 

w(ki) + w(k 2 ) = w(k 3 ) + w(k 4 ), 
ki + k 2 = k 3 + k 4 . 



(4) 



(5) 



Dynamical systems describing time evolution of the 
slowly changing amplitudes Aj of resonantly interact- 
ing modes can be obtained from (3), (4) or (3), (5) us- 
ing e.g. a multi-scale method. In a 3-wave system 
Aj = Aj(T), T = t/e and in a 4- wave system Aj = 

Aj(T), T = t/e 2 . The corresponding dynamical systems 
(in canonic variables) are written out below: 

iA l = ZA* 2 A 3 , iA 2 = ZA\A 3 , iA 3 = -ZA X A 2 ; (6) 



'iA x = VA* 2 A 3 A 4 + (wi - , 
iA 2 = VA\A 3 A 4 + (u) 2 - u 2 )A 2 , 
i A 3 = V*A\A l A 2 + (W3 - lu 3 )A 3 , 
i At = V*A 3 AiA 2 + (u)4 - W4) A 4 , 



(7) 



■<oi=EUi<yij\Aj\ 



2 ^S'il^l 2 ) ■ 



where the interaction coefficients Vu 



Vji = VII and 



V = V 34 2 are responsible for the nonlinear shifts of fre- 
quency and the energy exchange within a quartet corre- 
spondingly; (u>j — ujj) are Stokes-corrected frequencies. 
For very small nonlinearity, dynamical system (7) can 
be regarded in a simplified form, with u>j — ujj = 0, i.e. 
without nonlinear correction of frequencies. 

3-wave interactions dominate in a weakly nonlinear 
wave system if resonance conditions (4) have solutions 
and the coupling coefficients Z^0. Otherwise, the lead- 
ing nonlinear processes are 4-wave interactions. 

The following results hold likewise for resonances and 
quasi-resonances with small enough frequency mismatch. 



A. Intermittency, < e < 0.1 



II. D-MODEL 

Time evolution of a wave field in a weakly nonlinear 
wave system is described by a weakly nonlinear PDE of 
the form 

Lty) = ~eN{^) (3) 

where A is a nonlinear operator, < e <C 1 and L is 
an arbitrary linear dispersive operator, i.e. L(<p)=0 for 
Fourier harmonics ip — Aexpi[kx — w(k)i] with constant 



Excitation of a single mode in a 3-wave system gen- 
erates energy exchange within a resonance cluster only 
if this is the high-frequency mode w(k 3 ) from (4). In a 
4-wave system, excitation of a single mode generates en- 
ergy exchange only if it is the high-frequency mode w(k 3 ) 
in a Phillips quartet 

w(ki) +u;(k 2 ) - 2u;(k 3 ), k!+k 2 -2k 3 , (8) 

which is a special case of (5), [19]. Solutions of resonance 
conditions (4), (5) form a set of independent resonance 
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clusters. The form of a cluster uniquely defines its dy- 
namical system. 

Solutions of dynamical systems (6), (7) are known, 
[20, 21]; they describe periodic energy exchange within 
a resonant triad or quartet correspondingly. Resonance 
clusters of a more complicated structure may have a dy- 
namical system with periodic or chaotic evolution, [5] . 

In both 3- and 4-wave systems, resonant interactions 
are not local in fc-space; even more, in a 4-wave system 
with dispersion function u ~ k a , modes with arbitrary 
big difference in wavelengths can interact directly. In 
this case a parametric series of solutions of resonance 
conditions can be easily written out: 

'kf + k% = kg + k%, ki + k 2 = k 3 + k 4 , 
< ki = (k x , k v ), k 2 = (s,-k v ), (9) 
k,3 = (k x , ky), k 4 = (s, fcy), 

where s is an arbitrary real parameter (see Fig.l). 
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FIG. 1: Color online. Nonlocal interactions in a 4-wave sys- 
tem, u> ~ k a . Each couple of (red) dot-dashed lines with equal 
lengths correspond to specific choice of a parameter s. 

In any given 3- wave system, most of the modes are non- 
resonant. A non-resonant mode, being excited, does not 
change its energy at the slow time scale T. In the major- 
ity of 4-wave systems, each mode satisfies (5). However, 
excitation of a single mode does not generate resonance 
in the general case: the excited mode has to be the high- 
frequency mode in a Phillips quartet. 

B. D-cascade, e ~ 0.1 0.4 

D-cascade means a cascade computed in D-model by 
the ICEM method first presented in [6]. In both 3- and 
4-wave systems, D-cascades are generated by modulation 
instability (MI). Accordingly, the ICEM method can be 
applied for all PDEs in which MI has been established: 
focusing weakly nonlinear Schroedinger equation (NLS), 
[17]; modified NLS, [22, 23]; modified Korteweg-de Vries 
equation, [24, 25]; and Gardner equation, [26]. 

In both 3- and 4-wave systems, D-cascades are gen- 
erated by MI which is described as a particular case 
of the Phillips quartet (8) with uj\ = loq + Aw, W2 = 
luq - Au>, < Aw < 1: 

wi + co 2 = 2uj q , ki + k 2 = 2k . (10) 



The mode with frequency luq is called carrier mode. At 
each step of a discrete cascade, conditions (10) are satis- 
fied, with a new carrier mode generated from the previ- 
ous cascade step. 

Time evolution of the quartet (10) is studied in the 
frame of the nonlinear Schroedinger equation. The cor- 
responding time scale r = t/e 2 is called Benjamin- Feir 
time scale and is shorter than the time scale of resonant 
interactions. To understand this one has to take into 
account that small parameter e res < 0.1 yielding reso- 
nance interactions is in fact substantially smaller than 
Smi ~ 0.1 -r 0.4 corresponding to modulation instability: 
T = t/e 2 res > tje\ n = t. This fact is well established 
e.g. in the theory of wind generated oceanic waves, [29]. 

Conditions for MI to occur may be given as an insta- 
bility interval for initial real amplitude A and frequency 
ui of the carrier wave. For the NLS with dispersion rela- 
tion uj 2 — gk and small nonlinearity e ~ 0.1 to 0.25 the 
instability interval is described by 

< Auj/Akoj < y/2. (11) 

The most unstable mode in this interval satisfies the so- 
called maximum increment condition (in Bcnjamin-Feir 
form, [17]): 

Au;/ujAk=l. (12) 

For moderate nonlinearity, e ~ 0.25 to 0.4, the maximum 
increment condition reads (in Dysthe form, [22]), 

A(j/(uAk - ^ 2 A 2 k 2 ^j = 1. (13) 

Eqs. (12) and (13) each generate two chain equations 
(one for direct D-cascade and one for inverse D-cascade) 
describing the connection between the amplitudes of two 
neighboring modes in the D-cascade, under the following 
assumptions, [6]: 

(*) the fraction p of energy transported from one cas- 
cading mode to the next one depends only on the ex- 
citation parameters and not on the step number of the 
cascade; p is called cascade intensity; 

(**) modes forming a D-cascade have maximum insta- 
bility increment, i.e. a cascade is formed by the most 
unstable modes within the corresponding intervals of in- 
stability. This is a mathematical reformulation of the 
Phillips hypothesis that the spectral density is saturated 
at a level determined by wave breaking, [28]. 

In particular, (12) generates chain equations connect- 
ing mode n to mode n + 1 

uj n+ i = uj n + uj n A(cj n )k n , (14) 
u n+ i = u n - w n A(u n )k n , (15) 

for direct and inverse D-cascades correspondingly. This 
means that the D-cascades are formed by nonlinear fre- 
quencies depending on the amplitudes. 

From the chain equations various properties of D- 
cascades can be derived, including the form of the dis- 
crete and continuous energy spectra. 
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III. SURFACE WATER WAVES 

To demonstrate the wide range of the predictions 
which are given by our model we have chosen a classical 
example - surface water waves with dispersion function 
cj 2 = gk and small nonlinearity, e ~ 0.1 -j- 0.25. 

Before proceeding with our study we need to make an 
important remark on the terminology used below. Stan- 
dard vocabulary for discussing wave resonant interactions 
is "a 3- wave system" if (4), (6) are satisfied and "a 4- wave 
system" if (5), (7) are satisfied. Regarding resonance con- 
ditions for a Phillips quartet (8) one might formally con- 
clude that this is a system of three waves with frequencies 
wi, <x>2 and 2cj 3 . However, comparing the dynamical sys- 
tem for a 3- wave system (6) and the dynamical system for 
a Phillips quartet obtained from (7) by taking A3 = A4 
we can see immediately that these systems are different. 
Accordingly, a Phillips quartet may be referred to in the 
literature as a 4-wave system. 

In the text below we call the system (10) a 4-wave 
system though in the original papers whose results are 
interpreted using D-model this system is often called a 
3-wave system. Our terminology also allows us to avoid 
confusion while discussing cascade termination due to in- 
termittency in Sec. Ill C 3. 



A. Discrete and continuous energy spectra 

For determining D-cascade direction and scenarios of 
D-cascade termination we need first to compute the form 
of discrete energy spectrum. Detailed computation of D- 
spectra for various wave systems are given in [6] . For the 
convenience of reader below we outline this computation 
for surface water waves with small nonlinearity. 

All computations below are performed with chain 
equation (14) and yield energy spectra for direct cascade. 
Computations for inverse cascade should be conducted 
similarly but with chain equation (15); they are omitted 
here. 

Assumptions (*), (**) mean that E n = pE n _\ at any 
cascade step n , E n ~ A 2 n being the energy of the mode 
with amplitude A n . As the dispersion function in this 
case has the form ui 2 = gk, this allows to rewrite (14) as 



y/pA{u n ) = A(uj n + tJ 3 n A n /g) =J2^T ^l^ldY 

(16) 

(here notation A n = A(co n ) is used). 

Restricting ourselves to the first two terms of the Tay- 
lor expansion for the left-hand side of (16), we can obtain 
an ordinary differential equation and solve it analytically: 



ulX n A n /g + (1 - ^)A n = => (17) 



A n = g 



ui„ 



C,C = A -g 



w - 2 (18) 



Accordingly, the discrete energy spectrum for the direct 
cascade reads 



E n = E(co n ) ~A 2 n = g 2 



+ c 



1 2 



(19) 



where ui , A arc the excitation parameters and p = 
p(ui ,A ). 

The corresponding continuous energy spectrum E{ui) 
is computed as lim^oo \E n+1 - E n \/\ui n+1 ~ui n \ yielding 



E(tu)~2g 2 (1-VP) 



cai- 



rn 



In particular, the special choice of excitation parame- 
ters C = yields 



E(oj) ~ g 2 ui 5 



(21) 



which is the saturated Phillips' spectrum, [28]; this is also 
in accordance with the JONSWAP spectrum (an empir- 
ical relationship based on experimental oceanic data). 
Kinetic WTT predicts ~ in this case, [1, 2]. 



B. Cascade direction 

Combining chain equation and expression for the am- 
plitudes of the cascading modes we can study how cas- 
cade direction depends on the choice of excitation param- 
eters. 

For instance, for direct cascade oj n +i — co n > with 
C ^ the use of (14), (17), (18) yields 



< u n+ i -u) n = Lo n A(Lo n )/g 



(22) 



(Wg) ,-2 



c 



1 9 = (23) 



0--VP) 



u n + 



(l-yp) _ 2 

^0 - 9 7T — w 



"1/9 



(l-x/F) 



+ 



, (I -VP) -2 

A o - 9 ^ w 



"l/g 



A -g 



(i-VP) 



uJq > 



(24) 
(25) 
(26) 



9 (1 - VP) + [2A - g (1 - Vp)"o 2 \ "I > (27) 

As (1 — y/p) > 0, the range of frequencies forming direct 
cascade depends only on the sign of the expression 2A — 

5(1 - Vp^o 2 - 

An easy examination of (23), (27) shows how to choose 
excitation parameters A n , cj in order to observe direct 
cascade: 



if 2A >g(l- Vp)w ' 



(28) 



the only restriction on the range of frequencies forming 
direct cascade is trivial: uj n > w ; accordingly, only direct 
cascade will occur; 



if 2A <g(l- Vp)uj ; 



(29) 
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direct cascade will be observed for the range of frequen- 
cies Wq < w„ < w„ si where 



.9(1 -VP) 



g (1 - y/p)u - 2A 



(30) 



For simplifying further formulae we introduce here a 
small parameter e = ^o^o = Aou^/g and rewrite (30) 
as 



UJ i 



1 (i - VP) 
(l-y/p)-2e m 



(31) 



Physical meaning of the frequency u net is explained in 
Sec.IIIC2. 

Similar computations can be performed for inverse cas- 
cade, and also the case when both direct and inverse cas- 
cade are possible can be studied this way. In particular, 
for some choice of excitation parameters both direct and 
inverse cascade can be initiated simultaneously. This sce- 
nario is supported by wide range of experimental studies, 
e.g. [13-15]. 

All formulae (17), (18), (31) are given in terms of exci- 
tation parameters A$, cjo and cascade intensity p. This 
means that we should also compute p as a function of 
A , w , p = p(A D ,uj n ). This tedious computation will be 
given elsewhere. However, in the next section we give an 
example of the computation for a particular form of the 
solution (17). 

Notice that for studying predictions of the D-modcl 
in experimental data one can just measure ^Jp as the 
ratio of amplitudes of two consecutive cascading modes, 
Vp = A n+ i/A n , and apply formulae afterwards. 

C. Cascade termination 

1 . Breaking 

It was first shown in [30] that the amplitude of the car- 
rier wave may become so large that its steepness exceeds 
locally the maximum steepness of gravity waves yielding 
the onset of wave breaking. 

In order to demonstrate that this effect can be repro- 
duced in D-model, let us regard a particular solution of 
(17) with C = 0: 



(1- yp) _ 2 
A n =9 7T—U n \ 



As for this solution 



(l-yp) _ 2 

Ao = 9 oj = 

fp=(l-2 £o ) 2 

\A n =p n / 2 A a = (l-2eo) n A 0l 



(32) 

(33) 
(34) 



any choice of £o and A defines uniquely a cascade inten- 
sity p and the amplitude of the n-th cascading mode. 



It follows from (32), (33) that in this case all cascading 
modes have the same steepness e n = Sq, Vn: 



A h 



(1-Vp) 



£ . 



(35) 



This allows to compute the steepness e of the total wave 
packet at step n (before breaking) as 



(n+ l)e . 



(36) 



Accordingly, though the amplitudes of the cascading 
modes are decreasing, the steepness of the total packet 
is growing with an increasing number of cascade steps. 

For instance, direct computations demonstrate that if 
initial steepness e = 0.1, then after 3 cascade steps A ■ 
100%/^ s=s 0.5%. However, the total steepness of the 
wave packet is e = 4-0.1 ~ 0.4 and according to the 
Stokes criterion for the limiting steepness being about 
0.44, we conclude that mode A3 is about to break. A 
different choice of the initial steepness, say £0 = 0.05, 
yields the same total steepness e = 8 • 0.05 <~ 0.4 at 
the step n = 7 and cascading mode A-j contains about 
23% of the excitation energy while A ■ 100%/A 7 w 48%. 
Thus, varying excitation parameters one can predict the 
occurrence of breaking at the different cascade steps. 

Denoting limiting steepness of the wave package before 
breaking as £f, r , we conclude that the cascade terminates 
due to breaking if (nt r + l)£o = Sbr ~ 0.44, i.e. at the 
finite step n\, ri 



Ubr 



0.44/eo - 1. 



(37) 



At the end of this section we point out again that all 
results given by (33)-(37) are obtained for a specific form 
of solution of (17), namely, for C — 0. 

In the general case some results might be qual- 
itatively different: for instance, breaking may occur in 
the infinity rather than at some finite step. 

In this section we did not aim to present all possible 
formulae in their most general form but rather to demon- 
strate that growth of nonlinearity following by breaking - 
an experimentally well established phenomenon, [10, 13- 
15] - can be reproduced by the D-model. 



2. Stabilization 

If at some cascade step n s t the mode with frequency 
cj„ st is stable, then the condition (11) is not fulfilled, 
no additional mode can be generated and the D-cascade 
stops due to stabilization at some frequency Lo nat ■ 

From (11), (17), (18) it may be concluded that 







A nst U! nst k nst — 



(i-VP) 



(38) 
(39) 



Q= 0_^_ UJnst+Cuj Sj g ^ (40) 
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2 _ (i - Vp) 



/c = 



2 g(l-^p)iu- 2 -2A 

and for direct cascade stabilization occurs if 



; (1-^) 

(l-y/p)-2eo' 



(41) 



(42) 



which is in accordance with (31). 

It follows from (42) that direct cascade 

(a) stabilizes at the finite step uj < u nst < oo if 1 — 

Vp > 2e o; 

(b) stabilizes in infinity if 1 — ^/p = 2e ; then C = 
in (18) and corresponding continuous energy spectrum is 
Phillips spectrum <~ u;~ 5 (see Sec. IIIC 1); 

(c) stabilization does not occur if 1 — < 2eo while 
expression on the RHS of (42) becomes complex and has 
no physical meaning, i.e. stabilization conditions can 
never be fulfilled. 

Similar computations can be performed for inverse cas- 
cade. Though formally the termination conditions may 
allow the inverse cascade to be terminated at a negative 
frequency, this is physically irrelevant. This means that 
in a real physical system an inverse cascade terminates in 
some vicinity of zero frequency mode which might yield 
a substantial concentration of energy near zero frequency 
mode, also observed experimentally, e.g. [12]. 



3. FPU-like recurrence 

The fact that the long-time evolution of nonlinear wave 
trains of surface water waves may evolve in recurrent 
fashion (FPU-like recurrence), where the wave form re- 
turns periodically to its previous form, has been discov- 
ered experimentally and described in the pioneering pa- 
per of Lake et al., [16]. The next mile-stone step in the 
study of this effect has been performed by Tulin and 
Waseda in [13] where the authors refined the experimen- 
tal technique in a way that not only excitation frequency 
but also initial side bands and the strength of amplitude 
could be chosen. More experimental results can be found 
in [14, 15] and bibl. therein. 

In D-model, formation of a recurrent phenomenon (in- 
termittency) is due to formation of a cluster of resonant 
quartets, in the simplest case - an isolated Phillips quar- 
tet, (8). Its occurrence depends strongly on the form of 
the experimental tank. 

For some aspect ratio of the tank side lengths, inter- 
mittency can not occur as kinematic resonance conditions 



can not be satisfied. If for given aspect ratio, solutions 
of (5) exist, interaction coefficient 7^0 and initially ex- 
cited resonant mode(s) are modulationally stable, then a 
recurrence may be observed. 

Below we give a short list of experimental observations 
with their respective explanations: 

- no cascade is observed, rather recurrent patterns on 
the water surface are observed, [7]: 

initial steepness is too small to initiate modulation in- 
stability; 

- no intermittency is observed, rather a discrete cas- 
cade terminated by wave breaking, [10]: 

initial steepness is big enough to cause modulation in- 
stability and bJbr < w s t or stabilization is generally not 
possible for the chosen excitation parameters; 

- no intermittency is observed in the non-breaking 
regime, [31]: 

initial steepness is big enough to cause modulation 
instability, cascade terminates due to stabilization, i.e. 
uj s t < tobr, and the mode with frequency oj st is not a 
resonant mode in a resonant cluster possible for chosen 
experimental tank; 

- intermittency is observed in the non-breaking regime, 
[14, 15]: 

cascade stabilizes at the frequency uj st , the w st -mode 
is resonant mode and may excite a resonant cluster with 
another cascading mode. In particular, if u; st -mode and 
LUQ-raode form a resonance, complete FPU-like recurrence 
will be observed, [13-15]. If w st -mode forms a resonance 
with cascading mode with frequency uj =/= ujq, then par- 
tial recurrence will occur, with spectral peak being down- 
shifted to the frequency uj, [31]. 

- intermittency is observed at post-breaking stage, [13- 
15]: 

as essential part of the energy is lost due to breaking, 
amplitudes of newly excited modes may become modu- 
lationally stable and form a resonance with some of the 
previously excited cascading modes. This is only a quali- 
tative explanation, quantified prediction is an important 
separate topic which lies outside the scope of this paper. 
A possible theoretical scenario of the energy redistribu- 
tion at the post-breaking stage is developed in [13]. 

In this section we have shown how to use the chain 
equation to determine, depending on the excitation pa- 
rameters, the direction of the energy cascade and how 
the cascade will terminate. 

It should be noted that also the asymmetry of direct 
and inverse cascades as known from experiment, e.g. [13- 
16], may be deduced from the chain equation, [32]. 



IV. D-MODEL VERSUS KINETIC WTT 

During almost fifty years, kinetic WTT which requires 
a distributed initial state, was used to describe experi- 



ments using narrow frequency band excitation. This was 
considered legal, as the assumption was and still is that 
from the excitation frequency as a starting point quickly 
a distributed state will establish. The discrete part of 
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property 


D-model 


kinetic WTT 


assumptions 








1 


cascade origin 

in an S-wave system 


modulation instability, 
no dependence on S 


S-wave kin. eq., 
depends on S 


2 


initial state 


narrow frequency band 


distributed state 


3 


locality of interactions 


no assumptions 


necessary 


4 


existence of 
inertial interval 


no assumptions 


necessary 


5 


origin of 

cascade termination 


no assumptions 


dissipation 


6 


range of waves 
steepness 


< e ~ 0.1 -=-0.4 


< e < 0.1 


7 


cascade intensity 


is constant 


no assumptions 


8 


energy flux 


no assumptions 


is constant 


predictions 








1 


cascade is formed by 


nonlinear frequencies 


linear frequencies 


2 


spectrum form 

(a) 

(b) 


discrete and continuous, 
depends on 
the excitation 


continuous, 
does not depend 
on the excitation 


3 


transition from 
discrete to 
continuous spectrum 


included 


not included 


4 


direction of cascade 


included 


included 


5 


intermittency 


included 


not included 


6 


origin of 

cascade termination 


various scenarios: 

stabilization, 

breaking, 

FPU-like recurrence 


(see assumptions) 



TABLE I: Assumptions and predictions used in D-model and kinetic WTT 



the spectrum which was well observed in experiment was 
ignored in theoretical discussion, focusing on the contin- 
uous part of the spectrum. 

That this approach is not without problems was ac- 
knowledged within the community. As A. Newell noticed 
recently, "numerics seems to agree with the theory but 
experiments not", [33] (see also recent review [34]). In- 
deed, a distributed initial state as needed for applicability 
of kinetic WTT is easy to create in numerical simulations 
but not in laboratory experiments. 

Though D-model and kinetic WTT differ greatly in 
their assumptions and consequent range of applicability, 
sometimes the predicted form of continuous energy spec- 



trum is very close. To get more understanding which 
approach to apply in a given experimental setup we give 
the following comparison of the assumptions and predic- 
tions of D-model and kinetic WTT (short list is given in 
Table I). 

The crucial difference between descriptions of energy 
cascades in D-model and in kinetic WTT is the physical 
mechanism generating a cascade: modulation instability 
in arbitrary s-wave system versus s-wave interactions, 
S = 3,4,... . 

This means in particular that a D-cascade is generated 
by a mechanism which provides locality of interactions 
automatically. In the kinetic WTT the locality has to be 
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assumed, and no mechanism is suggested which allows to 
choose local interactions in wave systems where also non- 
local interactions are possible, as was shown in Sec. II A, 
Eqs.(9), and is also experimentally observed, [18]. The 
assumption of locality - only interactions among waves 
with close wavelengths are allowed - is basic in the kinetic 
WTT; without locality energy exchange among different 
scales k is possible and the energy spectrum can not be 
regarded MS (X function of only k. 

Another important point is that the influence of the 
excitation parameters on the form of the continuous en- 
ergy spectrum, observed experimentally, e.g. [12, 35-37], 
principally can not be included into kinetic WTT but is 
reproduced in D-model. 

One more considerable difference between D-model 
and kinetic WTT is the origin of cascade termination. 
In kinetic WTT this is always dissipation while in D- 
model various scenarios can be reproduced depending on 
the excitation parameters and direction of the cascade. 
D-cascades can terminate e.g. due to breaking, stabiliza- 
tion or formation of the Fermi-Pasta-Ulam-like recurrent 
phenomenon; all these effects are observed experimen- 
tally, [13-15]. 

Assumption (*) of D-model about constant cascade in- 
tensity, p = const, is absent in the kinetic WTT. This 
assumption is not substantial for D-model and can eas- 
ily be removed. Indeed, if cascade intensity at step n 
is p n 7^ const, chain equations (14), (15) do not change, 
while the ODE (17) and its solutions can be trivially 
rewritten by the changing p to p n . The only non-trivial 
change would be the construction of the transition from 
discrete to continuous energy spectra. Of course, the es- 
timates for determining cascade direction, termination, 
etc. should be recalculated and might get a more com- 
plicated form though not necessarily. For instance, all 
estimates made for the particular solution of (17) with 
C = remain valid while for so chosen excitation pa- 
rameters Aq,u>o cascade intensity is a constant defined 
by A ,w : 

C = => Ao-g Q—^- Uo 2 = => (43) 

p n = \Jl~ 2A u%/g = const . (44) 

Accordingly, transition from discrete to continuous spec- 
trum can be performed as above producing saturated 
Phillips spectrum. 

A wide range of experimental data shows that p = 
const in various wave systems and accordingly the dis- 
crete energy spectrum has exponential form, e.g. [38] 
and bibl. therein; this was our motivation for choosing 
constant cascade intensity in this presentation. 

Last but not the least. As it was shown in a recent 
experimental study of capillary waves, "from the mea- 
sured wavenumber-frequency spectrum it appears that 
the [linear] dispersion relation is only satisfied approxi- 
mately. (...) This disagrees with weak wave turbulence 
theory where exact satisfaction of the dispersion relation 



is pivotal. We find approximate algebraic frequency and 
wavenumber spectra but with exponents that are differ- 
ent from those predicted by weak wave turbulence the- 
ory", [39]. 

On the other hand, D-cascades are formed by the 
modes with nonlinear frequencies and not by the modes 
with linear frequencies as it is assumed in kinetic WTT. 

This is a manifestation of the very important differ- 
ence between cascades in the D-model and kinetic WTT. 
Cascades in the kinetic WTT are due to resonant inter- 
actions and therefore are possible at the time scales T or 
T with very small nonlinearity < e < 0.1. In D-model 
only intermittency is formed at these time scales while D- 
cascade occurs at the faster time scale r and for bigger 
nonlinearity £ ~ 0.1 -j- 0.4. 



V. CONCLUSIONS AND OPEN QUESTIONS 

In this paper we presented D-model which describes 
nonlinear wave systems with narrow frequency band ex- 
citation. It allows to reproduce in a single theoretical 
frame various nonlinear wave phenomena, in particular 
finite-size effects in resonators and formation of energy 
cascades. The cascades do not depend on shape or finite- 
ness of the interaction domain as they are triggered by 
the local mechanism of modulation instability. 

The main predictions of D-model can be stated as fol- 
lows: 

- Intermittency is formed by a set of distinct modes 
with linear frequencies; intermittency may occur in sys- 
tems with very small nonlinearity, < e < 0.1, at the 
slow time scales T or T; the underlying physical mecha- 
nism is resonant wave interaction. 

- An energy cascade is formed by a chain of distinct 
modes with nonlinear frequencies; a cascade may occur in 
systems with small to moderate nonlinearity, e ~ 0.1 -j- 
0.4, at the Benjamin-Feir time scale r; the underlying 
physical mechanism is modulation instability. 

- The discrete and continuous energy spectra of a cas- 
cade can be computed by the increment chain equation 
method, [6]; the form of spectra, cascade direction and 
scenario of cascade termination depend on the excitation 
parameters. 

- Various scenarios of energy cascade termination, 
known from laboratory experiments - stabilization, 
breaking and appearance of Fermi-Pasta-Ulam-like recur- 
rence, can be reproduced in D-model. 

As it was discussed in Sec. IV, all these predictions are 
quite different from those of kinetic WTT developed for 
wave systems with distributed initial state. In the lat- 
ter case an energy cascade occurs at the slow time scale 
of resonant interactions, is formed by linear frequencies, 
terminates (by assumption) always due to dissipation, 
etc. 

D-model explains known physical phenomena, as well 
as the results of individual laboratory experiments. In 
addition D-model makes predictions which may be eas- 
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ily verified in experiment, e.g. increasing the amplitude 
of excitation increases the distance between cascading 
modes in /c-space (direct consequence of the chain equa- 
tion), and others. 

It should be mentioned that within the wide range of 
excitation-dependent spectra predicted by D-model, the 
saturated Phillips spectrum w~ 5 has two special proper- 
ties. Firstly, as shown in Sec. Ill C 2, of all possible spec- 
tra only the Phillips spectrum does not stabilize after a 
finite number of cascade steps but in infinity (in fc-space) . 
Secondly, for the Phillips spectrum it is easy to prove that 
cascade intensity is constant (see (44)); for other spectra 
it is not known. What this physically means is presently 
under study. 

D-model may be refined in many ways, e.g. 

- in (17) just two terms of the Taylor expansion are 
taken to compute the energy spectrum; instead, one may 
regard the hierarchy of finite-order ODEs obtained by 
cutting off the Taylor expansion at 3, 4 and so on terms; 

- dissipation (depending on frequency) can be taken 
into account in the following way: cascade intensity p 
which describes the fraction of energy going from mode n 
to mode n+1 may be considered as a function increasing 
with frequency, p = p(w) ^ const. So stabilization of 
the cascade will occur earlier and also the form of the 
energy spectrum will change. 

Many more problems can be studied in the frame of 
D-model than have been mentioned in this paper. For 
instance, 

- Is it possible to use D-model for describing real-life 
phenomena where excitation parameters are not a priori 
known? 

Most naturally one might study the probability of var- 
ious initial states in a given situation and choose as input 
for the model either the most probable state or an aver- 
age state - for instance the known prevailing direction of 
the wind blowing over the ocean during a season. 

- Modulation instability plays a central role in the for- 
mation of extreme waves, e.g. [41-43]. Is it possible to 
use D-model to predict freak waves in the ocean? 

The Benjamin- Feir index (BFI), which is a ratio of the 



parameter of nonlinearity e to the relative spectral width, 
characterizes the evolution of an unidirectional wave field 
with a narrow spectrum. As either the frequency range 
or the directional spreading widen, the probability of ap- 
pearance of extremely steep waves decreases, [42, 43]. 
Using chain equation, one may e.g. to compute an upper 
estimate for BFI at each cascade step as a function of ex- 
citation parameters and to study characteristic behavior 
of this function. 

In the special case of surface water waves the Za- 
kharov equation is the model of choice. So it would be 
of great interest to compare the predictions of D-model 
with predictions of the Zakharov equation. 

Some results are known already, for instance, sideband 
asymmetry of Bcnjamin-Feir instability is established in 
numerical simulation with the Zakharov equation, [44]. 
Moreover, it was recently shown by M. Onorato, [45], 
that a D-cascade has a direct correspondence in the Za- 
kharov equation: the frequencies of cascading modes as 
determined in D-model form exact 4-wave resonances in 
the Zakharov equation with nonlinear Stokes corrected 
frequencies. 

This result is of the upmost importance as it opens a 
broad avenue for further studies of nonlinear wave sys- 
tems with higher degree of nonlinearity. The question 
is: 

Is it possible to compute energy cascades in nonlinear 
wave systems with distributed initial state using a new 
type of wave kinetic equation based on resonances of non- 
linear Stokes corrected frequencies, with bigger nonlinear- 
ity than is possible for applicability of kinetic WTT? 
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